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We present a scheme for creating quantum entangled 
atomic states through the coherent spin-exchange collision of 
a spinor Bose-Einstein condensate. The state generated pos- 
sesses macroscopic Einstein-Podolsky-Rosen correlation and 
the fluctuation in one of its quasi-spin components vanishes. 
We show that an elongated condensate with large aspect ratio 
is most suitable for creating such a state. 
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Quantum entanglement lies at the heart of the pro- 
found difference between quantum mechanics and clas- 
sical physics The entanglement between the states 
of space-like separated particles is the fundamental rea- 
son for the violation of Bell inequality, and causes many 
of the "paradoxes" of quantum physics. In recent years, 
there has been an interesting maturing of the discussions 
of entanglement away from the foundations of quantum 
mechanics and to "applications" in the emerging field of 
quantum information processing. 

A majority of the experimental realizations of quan- 
tum entanglement to date involve the creation of entan- 
gled photon pairs. Although ideal as carriers of quantum 
information, photons are however normally difficult to 
store for extended periods of time, in contrast to mas- 
sive particles. To overcome this difficulty, progresses 
have been made to generate correlated atom-photon pairs 
Recently, much attention has also been paid to 
quantum correlated atomic systems, particularly non- 
classical multi-atom states Q-Q, as these systems have 
important applications in quantum measurement beyond 
the "standard quantum limits" as well as in quantum 
computation. 

There have already been several proposals to create 
entangled atomic ensembles [|[^ and one of them has 
recently been demonstrated experimentally [0. All of 
these schemes rely on mapping the nonclassical proper- 
ties of electromagnetic waves, e.g., squeezed light, onto 
the state of an atomic system. In this Letter, we show 
that by taking advantage of coherent spin-exchange ul- 
tracold collisions, one can generate macroscopic atomic 
Einstein-Podolsky-Rosen (EPR) states from a spinor 
Bose-Einstein condensate without the need of nonclassi- 
cal light fields. 

We proceed by first giving the general idea of the pro- 
posed technique, and then turn to a more detailed the- 
oretical discussion. Our scheme is illustrated in Fig. ||. 
A spinor Bose-Einstein condensate consisting of a dilute 



F — 1 atomic sample is initially polarized such that only 
the spin-0 hyperfine ground state is populated at time 
t — Q. Binary spin-exchange interaction then convert the 
spin-0 atoms into pairs of spin-(±l) atoms. The irre- 
versibility of such a process is provided by shifting the 
energy of the spin-0 state above that of the spin-(±l) 
states (which can be achieved using the ac Stark shift 
provided by far off-resonant laser light [^). As a re- 
sult of this detuning, the phase-matching condition, i.e., 
conservation of momentum and energy, ensures that the 
resultant atoms in the pair move in opposite directions 
away from each other and escape the trap. Quantum en- 
tanglement results from our ignorance about which of the 
two escaping atoms is in the spin-(-l-l) state and which 
has spin-(— 1). 




FIG. 1. Entanglement scheme: A spin-0 condensate is ini- 
tially prepared. Spin-exchange interaction creates spin-(±l) 
atom pairs whose energy level is shifted below that of spin-0 
atom by an amount of fiS. This excess energy is transferred 
into the kinetic energies of spin-(±l) atoms which escape the 
trap. 

We now turn to a detailed analysis of this system. At 
t = 0, a condensate of iVo spin-0 atoms is confined in 
an optical dipole trap. An additional off-resonant opti- 
cal field is used to shift the energy of the spin-0 state 
above those of the spin-(±l) states by an amount fi5 (see 
Fig. 1^). The spatial wave function of the condensate, Lp{r) 
is determined by the stationary Gross-Pitaevskii equa- 
tion. 

At i > 0, the spin-(±l) states start being populated 
by the spin-exchange interaction 

H = Kj dr^l^{r,t)ipl^{r,t)iPo{r,t)Mr,t) + h.c. (1) 

where Xa is a constant related to the s-wave scattering 
lengths associated with the hyperfine levels involved, and 
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^pa is the boson annihilation operator for spin-a atoms. 
The effect of atomic recoil during this process is to trans- 
fer the excess energy hS into the kinetic energy of the 
spin-(±l) atoms. Therefore, for the short time scale 
where the propagation of (±1) atoms can be neglected, 
we may expand the boson field operators as 



(2) 
(3) 



where ujg = ?i|qp/ (2m), and the operators {c^} obey the 
boson commutation relations [c^, cj,] = ^^.i,. With these 
expansions, Hamiltonian (fil) may be reexpressed as 



(4) 

where Aq^q, = {ujq + ujq, - 2(5), K = \aV'^/{2'Kf with V 
being the quantization volume, and 



P(q,q')= j dr|(^(r)|'*e-*('i+i')-''. 



(5) 



Eq. (Q) is reminiscent of the Hamiltonian describing para- 
metric down conversion processes in nonlinear and quan- 
tum optics. As is well known, these processes lead to 
squeezing and to the generation of entangled photon 
pairs. 

For short enough interaction times, the population of 
the sidemodes (±1) remain small compared to Nq. In this 
regime, we neglect the depletion of the spin-0 state and 
treat cq as a c- number cq such that |coP = A^o- We can 
furthermore neglect those terms in the Hamiltonian (^) 
that describe atom-atom interactions involving only the 
spin-(±l) states. Under these assumptions, the Heisen- 
berg dynamics of the operators c±i^q resulting from the 
interaction Hamiltonian M) simplifies to 



c+i,q = -iK j rfq' p(q, q') e 

^c_i,q = -iK/ fiqV(q,q')e 



--l,q' ^0' 



'^+l,q'C0- 



(6) 
(7) 



where we have defined the gain parameter 

G^^2nK^J dq'|p(q,q')p5(A,,,0 (10) 

and the noise operator 

fl{t)^^iK4 J dqV(q,q')e'^-'*cLi^q,(0), 

whose correlation functions are given in the Markov ap- 
proximation by 

{flit)Ut')) = o, 

(/q(i)4'(0)=A^oGq<5(q-q')<5(i-0- 

It is this noise operator that triggers the populat- 
ing of spin-(-l-l) state from quantum fluctuations. In 
deriving Eq. (|9|), we have used the approximation 

p(q,q')P*(q',q") ~ l/o(q,q')P^(q - q") and neglected 
the principal part associated with the definition of the 
(5-function. 

Following a similar procedure, we can derive the equa- 
tion of motion for c_i.q as 



^C_l^q=^GqC_l,q-f5t(t), 



(11) 



where 



To solve these equations, we first formally integrate 
Eq. (0) to get 

c-i,q(t) = c_i,q(0) - iKcg y dq'p(q,q')y dre^^-^'" 

« c_i,q(0) - iKcg J dqV(q,q')^(A,,,')fili,q(0 (8) 

where the Markov approximation has been invoked. In- 
serting Eq. (H) into Eq. (^), we obtain 



|c+l.q=^GqC+l,q + /t(t) 



(9) 



aiit) = / dqV(q,q')e*^-'*c^i,q,(0)c2. (12) 

At this level of approximation, which neglects as we 
recall the depletion of the spin-0 mode, the Heisenberg 
equations of motion (||) and (^ij) are linear. They can 
readily be integrated to give 



2+l,q(i)-aq(i)c+l,q(0)+ / dr (t ) /t (t - T ) (13) 

Jo 

5-l,q(t)=aq(t)c_l,q(0)+ / dr ^q (t) .9 1 _ r ) (14) 

Jo 

where Gq{t) = exp(A"Q Gqt/2). From these we can calcu- 
late the population in modes {±1, q}: 

A^±l,q = (4l.qC±l,q) = eXp{N^Gqt) - 1 

It is also straightforward to calculate the correlation func- 
tion 

Cq,q' = (C-l,qC+l^q/} 

= -»5,WcSp,(c,i|g2_£!!i£!. ,15, 

The fact that the ±1 modes are correlated implies that 
the two spin states (±1) are entangled. It is obviously 
desirable that a spin-(— 1) atoms with momentum hq be 
correlated to a spin-(-l-l) atom with well defined momen- 
tum hq' . 
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From the definition (||) of p(q, q'), we conclude that 
as long as the spatial size of the condensate wave func- 
tion is much larger than the reciprocal length l/|q| and 
l/|q'|, p(q, q') is approximately proportional to a delta- 
function, p(q, q') — > (5(q -t- q')- In other words, under 
this condition the two correlated atoms resulting from a 
spin-changing collision move in opposite directions. Ad- 
ditionally, the particle momenta |q| and |q'| have to sat- 
isfy the conservation of energy condition ujq+ujq'~2S w 0, 
a condition that can be met for a large light shift hS. We 
observe that in addition, a large energy shift h6 is also re- 
quired to produce spin-(±l) atoms with sufficiently large 
kinetic energy to escape the trap. As already mentioned, 
this is required to prevents them from undergoing a col- 
lision resulting in a pair of spin-0 atoms. 
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FIG. 2. The gain parameter Gq versus 6, for a Gaussian 
and cylindrically symmetric condensate wave function of the 
form (p{r) oc exp[-z^/(2crf ) - {x^ + y^)/(2cri)]. 9 is the 
angle betwee n q and the z-axis. In the calculation, we set 
9 = IqI = y^2TnS/h. Curve 1; — 10, q — 10; curve 2; 
CTz = 10, g = 20; curve 3: = 10, q — 40; curve 4: = 20, 
q = 40. The units for and q are a± and l/o-±, respectively. 



In general, it is not sufhcient to just produce entan- 
gled atomic pairs. Rather, one would like to subsequently 
store them, e.g. in a dipole trap. It is desirable for this 
purpose to achieve a high degree of directionality in the 
generated atoms, so that they have a narrow enough an- 
gular distribution. To see how this can be achieved, let 
us take a closer look at the gain parameter Gq appear- 
ing in Eq. (|lO|). The expression of Gq is reminiscent of a 
similar gain parameter encountered in the study of super- 
radiant scattering from a condensate [|l3[|l4| ] . It has been 
shown in that context that for a spatially anisotropic 
condensate, the largest gain occurs along the longest di- 
mension of the condensate . The same conclusion can 
be reached in the present case. Fig. || illustrates the gain 
along different directions for the case of a cylindrically 
symmetric condensate, for various light shifts hS and as- 
pect ratios. For simplicity, we choose q — — y^2m5/h, 
and assume that the condensate has a Gaussian shape. 
Fig. 2 illustrates quite clearly that a smaller angular 
distribution of emitted atoms is obtained for larger as- 



pect ratios and larger q. Thus for a strongly elongated 
cigar-shaped condensate, the matter-wave modes along 
the long axis, which have the largest gain coefficient Gq, 
will typically deplete all the condensate atoms before the 
population of the off-axis modes can significantly build 
up. As a consequence of mode competition, the emis- 
sion of the spin-(±l) atoms is therefore largely confined 
to two narrow cones at the two ends of the cigar-shaped 
condensate. 

From this discussion we conclude that in order to ex- 
perimentally realize the proposed scheme, one should first 
create an elongated spin-0 condensate with a large light 
shift hS. Spin-exchange interactions then generate pairs 
of spin-(±l) atoms fiying in opposite directions along its 
long axis. These atoms can be subsequently captured by 
two traps located at opposite sides of the original trap. 
Eventually, the spin-0 condensate is depleted, with two 
new ensembles of pair-wise entangled atoms stored in- 
side the side traps. We emphasize that although each 
trap contains both spin-(-l-l) and spin-(— 1) atoms, these 
cannot undergo subsequent spin-exchange collisions to 
produce spin-0 atoms, since this process does not satisfy 
momentum- energy conservation. 

The spin-(±l) atoms being created in pairs, we know 
for sure that taken together, the two ensembles must con- 
tain equal number of spin-(-l-l) and spin-(— 1) atoms — 
although how many spin-(-l-l) and spin-(— 1) atoms are 
in each ensemble is unknown. In the Schrodinger picture, 
such a state may be expressed as 



N/2 
m=-N/2 



N 



N 



(16) 



where N/2 is the total number of atoms in each of the two 
"left" and "right" side traps, labeled by / and r, respec- 
tively. The integers m and —m represent the difference 
in numbers of atoms in the spin states (+1) and (—1) in 
each of the two traps. 

Introducing the z-component of the quasi-spin opera- 
tor 



W = N 



+1 



N 



(i) 
1' 



where is the number operator for state-(±l) in en- 
semble i and i — l,r, we have that 



N 



N 



Since the explicit expressions of the coefficients Um in ( p^ ) 
are unknown, so are the expectation value and variance 



for L 



However, a simple calculation shows that 

(4) - 0, 
(AL,)2 = 0, 



where Lz 



HI) 



Lz is the ^-component of the total 



quasi-spin operator. Hence, although the variance of the 
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may be large for the individual ensembles, the vari- 
ance for the whole system vanishes. In other words, taken 
as a whole the two ensembles represent a maximally spin- 
squeezed state. This should be contrasted to the case of 

independent atoms in the state (| + 1) + | — 1))^, foi' 
which one finds (AL^)^ = N/4. 

We note that if we randomly pick one atom each from 
the two side-traps for an atomic ensemble prepared in 
state (Iq), then their degree of entanglement is only of 
order iJN. This is because although the atoms are cre- 
ated in pairs, we cannot tell which particular pairs of 
atoms are entangled. It is only through the collective 
spin measurement that the quantum entanglement can 
be revealed. The observation of such an macroscopic en- 
tanglement can be carried out with the technique of spec- 
troscopic detection of collective spin noise at the quan- 
tum level described in Refs. |lO|Jl^ ]. In practice, the state 
( p^ ) has to be averaged over the statistical distribution 
of the total particle number N. However, as noted in 
Ref. 1^ , such fluctuations do not affect the entanglement 
significantly for large numbers of atoms. 

In conclusion, we have proposed and analyzed a simple 
scheme to create a macroscopic EPR-correlated atomic 
state. Such a state possesses a nonlocal entanglement 
and is maximally squeezed, in the sense that the fluctua- 
tions of the z-component of its quasi-spin vanish. Hence 
we believe that this system will have important applica- 
tions in precision measurement as well as in fundamental 
physics such as the test of nonlocality in macroscopic 
quantum systems. Our study shows that an elongated 
spinor condensate with large aspect ratio and large en- 
ergy difference between spin-0 and spin-(±l) states is the 
best candidate to create such a state. The correlations 
between the atomic ensembles arise from the nonlinear 
atom-atom interaction amongst the condensate atoms. 
This distinguishes our work from other proposals with a 
similar goal, where the correlations between atoms are 
transferred from EPR-correlated light fields. As a con- 
sequence, our scheme can deal with strictly ground state 
hyperfine atomic states. This is of considerable advan- 
tage, since the entanglement of the kind described here is 
therefore robust against decoherence and immune from 
the quantum fluctuations caused by the electromagnetic 
vacuum field modes, which limit the degree of entangle- 
ment and spin squeezing 

Note: Upon completion of our work, we noticed a 
preprint paper by S0rensen et al. | [T^ ] in which the possi- 
bility of creating squeezed spin state with Bose conden- 
sates is investigated. Our work differs from theirs as the 
atomic ensembles we studied are spatially separated with 
nonlocal EPR correlation, while theirs does not possess 
this property. 
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